Quantum theory of Synchronously Pumped type I Optical Parametric Oscillators: 

characterization of the squeezed supermodes 
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Quantum models for synchronously pumped type I optical parametric oscillators (SPOPO) are 
presented. The study of the dynamics of SPOPOs, which typically involves millions of coupled 
signal longitudinal modes, is significantly simplified when one considers the "supermodes", which are 
independent linear superpositions of all the signal modes diagonalizing the parametric interaction. In 
terms of these supermodes the SPOPO dynamics becomes that of about a hundred of independent, 
single mode degenerate OPOs, each of them being a squeezer. One derives a general expression 
for the squeezing spectrum measured in a balanced homodyne detection experiment, valid for any 
temporal shape of the local oscillator. Realistic cases are then studied using both analytical and 
numerical methods: the oscillation threshold is derived, and the spectral and temporal shapes of 
the squeezed supermodes are characterized. 
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I. INTRODUCTION 

Mode-locked trains of pulses, or frequency combs, have 
at the same time the coherence properties of c.w. lasers 
and the high peak powers of pulsed lasers. They are po- 
tentially perfect tools for generating non-classical states 
of light, as they are are at the same time high quality light 
sources, free of excess noise, and intense sources able to 
induce strong nonlinear effects, and therefore to generate 
strongly non-classical states of light, such as squeezed or 
quadrature-entangled states. 

Frequency combs have been used in many quantum op- 
tics experiments, and have efficiently produced non clas- 
sical light, either in X (2) 0,0,1 or X (3) @, 0, S, 0, S S 
ITol . ITlj | media, but so far in a single-pass configuration in 
the nonlinear medium. In this configuration, one needs 
very high peak powers, and the system loses somehow 
its potential high quality in terms of pulse to pulse co- 
herence and transverse profile. Mode locked lasers have 
also been used to efficiently generate squeezed states in 
optical fibers [H, [H[ . The system has the drawback of 
generating non- minimal states states of light, because of 
the excess noise due to Brillouin scattering in the fiber. 

We have recently proposed [lij to use synchronous 
optical cavities to recirculate the light in the nonlin- 
ear medium, thus enhancing further the nonlinear effects 
and imposing the cavity mode structure to the generated 
non-classical field. This can be done by building "syn- 
chronously pumped OPOs" or SPOPOs. In a SPOPO 
the cavity round-trip time is equal to the delay between 
successive pulses of the pumping mode-locked laser, so 
that the effect of the successive intense pump pulses 
add coherently, thus reducing considerably its oscillation 
threshold. In [lj], a large squeezing effect was predicted 
in some "supermodes" , which are well defined linear com- 
binations of signal modes of different frequency, but not 
studied in detail. In the time domain these supermodes 
correspond to trains of pulses of different waveforms, or- 
thogonal each other. The purpose of this paper is to 



precise the quantum model used to predict the effects 
and to investigate in a detailed way, through analytical 
or numerical methods, the potentialities of the system in 
realistic situations 

SPOPOs have already been implemented as efficient 
sources of tunable ultra-short pulses [H, [I?], [H, 111, 
I20L HH their temporal properties have been theoretically 
investigated [22|, [H, [24|, and actively mode-locking of 
OPOs has been recently achieved [25[. 

The decomposition of a pulsed field in terms of a ba- 
sis of normal modes, similar to the supermodes we con- 
sider in this paper, has been introduced in different con- 
texts for a complete quantum characterization of cither 
the pulsed squeezed light generated by parametric down 
conversion [2y, [2?J or solitons in optical fibers 28] . Such 
approaches are strongly connected with the Schmidt de- 
composition of two-photon states for the characterization 
of pairwise entanglement [29|, [3(| and the Bloch-Messiah 
reduction of any optical circuit characterized by a lin- 
ear input-output relation (3l| . In this context Menicucci 
et aJ. [32|] proposed optical frequency combs as scalable 
resources for quantum computation. 

The article is organized as follows: we present first 
the model that we will use. The system turns out to be 
characterized by a real and symmetric matrix C, which 
contains all the information about the effective nonlinear 
interaction. The eigensystem of C is thus of special rel- 
evance and is studied in Section III, where the SPOPO 
threshold and several general properties of the spectrum 
of C are addressed. An analytical approximation to the 
diagonalization is also given that allows a better insight 
into the general trends as parameters are varied. In Sec- 
tion IV it is shown that the introduced eigenmodes or 
supermodes are squeezed, the corresponding eigenvalues 
determining the amount of squeezing, which can be mea- 
sured in a balanced homodyne detection experiment that 
uses as the local oscillator (LO) a field with the same 
spectrum as the desired supermode. In Section V one 
then studies the squeezing properties of SPOPOs in two 
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realistic cases, corresponding to BIBO and KNDO3 crys- 
tals, using an appropriate scaling property of the diago- 
nalization problem. Finally, an appendix details the case 
of the singly resonant SPOPO. 



II. THE SPOPO MODEL 

A. Evolution equations for the operators 

We consider quasi-degenerate collinear type I interac- 
tion, by means of which the pumping frequency comb, 
at frequencies around 2ujq, is converted by a nonlinear 
•^(2) crystal into multimode signal radiation at frequen- 
cies around u/q, and vice- versa, where 2loq and ujq are the 
two frequencies at which perfect phase matching occurs. 
This implies that one has n (2u>q) — n (ojq) = uq, where 
n is the crystal refractive index. The nonlinear crystal 
is placed inside a high finesse optical cavity of length L, 
which is assumed to be dispersion compensated by in- 
tracavity dispersive elements, so that all cavity modes 
around the frequency loq are equally spaced by a com- 
mon free spectral range f2, which is made equal to that of 
the pumping laser, thus warranting the synchronization 
of the pump to the OPO cavity. This ensures that the 
pulse-to-pulse delay of the pump beam coincides with the 
cavity round-trip time and successive pump and signal 
pulses superpose in time, thus maximizing the strength 
of the interaction. Hence the external pump mean field, 
which is a phase-locked multimode coherent field, can be 
written as 



which extends from z — —lj2 to z = +1/2, one can write 



-Eext (t) 



2enc ^ 



(1) 



P is the average laser irradiance (power per unit area), 
a m is the normalized (X) m \ a m\ — 1) complex spectral 
component of longitudinal mode labeled by the integer 
index to, and to = corresponds to the phase-matched 
mode. As we will be concerned with femtosecond lasers 
with pulse durations around lOOfs , the number of pump 
modes will be typically on the order of 10 4 — 10 5 . 

Two possibilities for pumping can be used: either (i) 
the pump also resonates inside the cavity (doubly res- 
onant case), which requires in addition dispersion com- 
pensation at the pump spectral region, or (ii) the cavity 
is transparent for the pump (singly resonant case), a case 
which is free from the previous restriction and thus more 
amenable for experimentation, at the expense of a higher 
threshold, as we will see. We detail more the latter case 
in the appendix at the end of this paper. We will limit 
here our analysis to non-chirped pumps as chirping re- 
quires a more general treatment which will be presented 
elsewhere. 

As the finesse of the cavity is assumed to be high, the 
intracavity signal field operator E s can be written as a 
superposition of cavity modes. Inside the crystal, 



E s (z, t) = ^2 i£s,mSm(t)u m (2) ' 



H.c. 



(2) 



where 



ujq + mf2, s m (t) is the annihilation oper- 



ator for the TO-th signal mode in the interaction picture, 
verifying standard boson commutation relations 



(*)>4 (*)] = <Wj 



(3) 



u m (z) is the spatial profile of mode m, equal to e lk "' mZ 
in the case of ring cavities, while for linear cavities it is 
equal to sin [k s ^ m (z + L/2)], where 



k s ,m — k (oj s ^ m ) 



n(u) Stm )u s 



(4) 



is the corresponding wavenumber. Finally 8 SjTn is the 
single photon field amplitude, whose value depends on 
the type of cavity. For a ring cavity 



£(ring) 



2e n (u s , m )A s L' 



(5) 



where ^4 S is the transverse area of the signal field, while 
for a linear cavity 



c(linear) _ /oc(ring) 



(G) 



Note that we are writing the field as a superposition of 
plane waves, but the treatment is approximately valid for 
Gaussian beams provided that the crystal is placed at the 
beam waist and the Rayleigh length is much longer than 
the crystal length I. In this case ^4 S = irw 2 with uu s the 
beam radius. Similarly we have A p — ttw 2 for the pump 
transverse mode. 

The interaction Hamiltonian Hi describing the para- 
metric interaction in the nonlinear crystal is given as 
usual by: 



+1/2 

Hi = -Ai f dz 

-1/2 



E p (z,t)P p (z,t) + E s (z,t)P s (z,t) 



(7) 

where P s (z, t) and P p (z, t) are the nonlinear electric po- 
larization at signal and pump frequencies, and Ai ac- 
counts for the effective area of interaction correspond- 
ing to the three-mode overlapping integral across the 
transverse plane and, for Gaussian beams, it is given by 
A^ 1 = A^ 1 + 2A~ 1 . The calculation of the Hamiltonian 
depends on the type of configuration (singly- or doubly 
resonant). Here we consider the simpler case of a dou- 
bly resonant SPOPO and leave the details of the singly 
resonant case for the Appendix. 



1. Doubly resonant SPOPO 

In this case an expression for the intracavity pump field 
operator E p analogous to @, now centered around 2ujq, 
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can be used and the following expression for H\ in the 
rotating wave approximation is obtained: 



Hi = 2ie xlAi ^£ s , m £ s ,g£p 



X 4W 4 {t)Pm+q (t)+H.C, 



(8) 



where x is the relevant nonlinear susceptibility, p m (t) 
and p' m (t) are pump boson operators (to = denotes the 
phase matched mode) verifying [p m (t) , p^ (t)] = 5 m ,„, 
and £ Pim is as £ s , m with the substitutions u) 8 ,m —> Up,m = 
2luq + mQ and ^4 S — > A p . The phase-mismatching factor 
of the crystal, f m ,q, is given by: 



fm,q — 



u m,q 



IJ m,q 



b m ,q being the phase-mismatch angle 

<f>m : q = — 



(^p,m+g ^s,m ^s,<j) ^- 



(9) 



(10) 



Making use of the standard input-output formalism 
of optical cavities the following set of Heisenberg equa- 
tions for the signal annihilation operators s m is derived 
straightforwardly: 



K ^ ] fm,qSqPm+qi (H) 



where the cavity damping rate %, or cavity linewidth, is 
equal to T s <C 1 is the transmission factor of the 
single cavity mirror at which losses are assumed to be 
concentrated, the coupling constant k is given by 



Qxi- 



f~AZ \n Lj 



(12) 



and Q is a factor depending on the geometry of the cavity 
that amounts to 2 for the ring cavity and to \[2 for the 
linear cavity. 

Analogously, the evolution equations for the pump an- 
nihilation operators p q are: 



dp q 
dt 



I rn^q^rn^q—m-! 



(13) 



where 7 P is the cavity damping coefficient evaluated at 
pump frequencies. 

To get these simple equations, we have assumed that 
£s,m — £s,o Vm, and neglected the dispersion of the non- 
linear susceptibility, which is a very good approximation 
as far as the pulses bandwidth is not too large [35J. The 
"in" operators correspond to quantum fields entering the 
cavity through the coupling mirror. We consider the 
case where the input signal field is the vacuum, and the 
input pump field a coherent state. We have therefore 
(s in ,m (*)) = 0, (p in ,g (f)) = p cxt , q , and the following cor- 
relations 

(An,m(t),pL, ro <(*0) = (Sln,m(*),4,m'(0 

= 5 m , m >S(t--e), (14) 



with the notation (a,b) = ({a — (a)) — \J^j\ the 
rest of correlations being null. The mean input field p cx t,q 
is related to the a q and P coefficients introduced in ([T|) 
by: 



Pcxt,q 



n A v P 
— ( 

2Hujq 



(15) 



2. Singly resonant SPOPO 



As detailed in the Appendix, in the case where only the 
signal field is resonating into the cavity the evolution of 
the signal field annihilation operators is given by an ex- 
pression identical to (|11[) , but now the pump annihilation 
operators are given by 

p m (t) =gJ^(til(t)+P^lM) + 

C ' 

n 

We see that these operators contain two contributions. 
The first one corresponds to the free-field part that is 
described by means of two independent boson operators 
Pin m W associated to the fields impinging the cavity from 
both the directions labeled with the superscripts (±). We 
consider the unidirectional pumping case where the input 
pump field propagating from left to right (labeled with 

(+)) is a coherent state with a mean value <rf <?£>(*)) = 
Pcxt,q, while the input pump field propagating form right 
to left (and labeled with (— )) is the vacuum (j>i~ q (t)) = 0. 
The mean input field, p ex t,q, is still given by Eq. (|15p . 
and the only non-null correlations are: 



j m, rn 



S m ,n S(t — t') . 



(17) 



B. The SPOPO below threshold 



Below threshold signal modes have a zero mean value, 
whereas the pump field is characterized by a huge am- 
plitude. One can therefore use a linearization procedure 
for the quantum fluctuations, which amounts to setting 



(P 



m+q/ 



In the doubly resonant case 
^/2/^/ p p ext<m+q , as given 



Pm-\-q 

(Pm+q) = y/Vlp (Pin,rn+g 

by (|13p . while in the singly resonant case (p m +q 
G \fW c (Piljm+q) = G-y/L/cpe^t,m+q, as given by CE 



ext,m-t-g. 

In both cases the final equations for the signal field ani- 
hilation operators are identical: 



~dt 
where 



-IsSr. 



\fPjPo, 



(18) 



(19) 
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is a pump amplitude parameter, and Po is an important 
scaling parameter for the pump, which can be shown to 
be the threshold value for the pump in the c.w. regime 
(single mode pump configuration). In the optimized con- 
figuration for the pump focussing (A p = A s /2), it is equal 
to: 

e c 3 n 2 T 2 



p = n - 



(20) 



2 (xM 

where IIo is given in Table [I] for the singly or doubly res- 



no 


Doubly resonant 


Singly resonant 


Linear cavity 


Tp/16 


1 


Ring cavity 


Tp/4 


1 



TABLE I: Coefficient for retrieving the cw threshold Po in the 
different experimental situations considered. 

onant configurations and ring or linear geometries. The 
differences arise from the fact that the doubly resonant 
case presents an intracavity pump power enhancement 
factor of 4/Tp with respect to the singly resonant case 
and the crystal is used twice in a linear cavity as com- 

(doubly) / p(singly) 



pared to the ring one. Hence the ratio P ( 



o 



IK 



equals the (very small) pump transmission factor of the 
doubly resonant cavity. By way of example, if we con- 
sider a singly resonant SPOPO (the geometry does not 
matter) based on a BIBO crystal [33| with a thickness 
of I = 100 /im and pumped at 0.4 /im, we obtain refer- 
ence irradiances Po of approximate values 14, 344, and 
1400 MWcm- 2 for T s = 0.01, 0.05, and 0.1, respectively. 
For a typical pump beam radius of 70 /im these irradi- 
ances lead to pump powers equal to 2, 53, and 212 kW, 
respectively. 

The key point for the following analysis is the fact that, 
in equations (fT8|) , the parametric coupling between the 
different signal modes is linear. It is characterized by a 
matrix C, with matrix elements: 



sm ( 



q^m+q 



m,q 



im+q • 



(21) 



When necessary, the phase mismatch angle 4> m .q, Eq. 
(|10[) can be computed using a Taylor expansion around 
2wo for the pump wave vectors fc Pjm and around loq for 
the signal wave vectors k SiTn , 

4> m>q ~f3 1 {m + q) + f3 2p (to + qf - f3 2s (to 2 + q 2 ) , (22) 
where 



1 



ft = -n (fcl - k' s ) I, 



2" v p 
1 

_( 

4" 

(3 2s = -tfKi, 



(3 2p = -tfk' p % 



(23) 
(24) 
(25) 



are dispersion coefficients, and k' and k" are the first 
and second derivatives of the wave vector with respect to 
frequency. Note that the matrix C depends on the cavity 
characteristics only through the free spectral range f2. 



III. SPOPO DYNAMICS FOR THE MEAN 
FIELDS. DETERMINATION OF THE SPOPO 
THRESHOLD 

Before calculating the quantum fluctuations of the 
SPOPO we analyze first the dynamics of the mean values 
of the operators, which is obtained by removing in Eq. 
(|11[) the input noise terms and replacing the operators 
by complex numbers: 



-IsSr. 



(26) 



The solution to Eqs. (|26[) is of the form 

s m {t) = S Km e Xkt , (27) 

where k is an index labelling the different solutions, and 
the parameters Sk, m and Xk obey the following eigenvalue 
equation: 

Afc5fe, TO = —J S Sk,m + 7sO" £ C-m,qSh,q- (28) 

Q 

As matrix C is both self-adjoint and real, its eigenval- 
ues Afc and eigenvectors Lk, of components Lk m defined 
by 



(29) 



are all real. As 7 S and a are also real, it is evident that 
two sets of solutions to Eqs. ([28]) exist, namely Sj^H = 
Lk t m and ^ = iLk,mi with corresponding eigenvalues: 



A 



(±) 



7, (-1 ± o-Ak) • 



(30) 



Let us label by index k = the solution of maximum 
value of | Afc | . When a |Ao| < 1, all the rates are neg- 
ative, which implies that the null solution for the steady 
state signal field is stable. For simplicity of notation, we 
will take A positive in the following, which is a com- 
mon situation as shown below (361 ■ Hence An + ^ is the 



largest eigenvalue and A 



sets the SPOPO oscilla- 



tion threshold, which then occurs when the pump param- 
eter a takes the value 1/A , i.e. for a pump irradiance 
P = Pthr equal to: 



Pthr = Po/A 2 ,. 



(31) 



The exact value of Ao, and therefore of the SPOPO 
threshold, depends on the exact shape of the phase 
matching curve and on the exact spectrum of the pump 
laser. As will be shown in Sec. V, the theoretical SPOPO 
threshold can be extremely low, of the order of the cw 
single mode threshold divided by the number of pump 
modes. 
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Let us now define the normalized amplitude pumping 
rate r by 



thrj 



or r = crA , so that the threshold occurs at r 
eigenvalues X k become 



A 



(±) 



7s 



-1 ± r- 



A 



(32) 
1. The 

(33) 



We will call supermodes the set of Sk, m values for a 
given k, which corresponds physically to the different 
spectral components of the signal field, and critical su- 



permode Sj^ , 



the one associated with A 



(+) 



which 



is the eigenvalue changing its sign at threshold. Above 
threshold, this critical mode will be the "lasing" one, i.e. 
the one having a non-zero mean amplitude when r > 1. 
Note that the supermodes are independent of pump (Eq. 
,but not the eigenvalues (Eq. |3"0)) . 
We note that the supermode in quadrature with re- 



spect to the critical one, S t 



(-) 
o 



iS, 



(+) 
o 



has an associ- 



ated eigenvalue Aq = —2j s at threshold, Eq. (|33|) with 
r = 1 , which is the lowest eigenvalue below- or at thresh- 
old. This property is obvious: Should ArT^ < — 2^/ s for 
some k, then r^j- should be larger than 1, what is incom- 



patible with the fact that 



< 1 by definition and the 



condition r < 1 . The fact that there exists an eigenvector 
whose damping rate (Aq ' = — 2j s at threshold) is twice 
that of the passive cavity has important consequences on 
the squeezing properties of the SPOPO, as it occurs in 
other OPO configurations [53] • 



IV. QUANTUM FLUCTUATIONS OF THE 
SPOPO BELOW THRESHOLD 



A. Fluctuation spectrum for the supermodes 

We can now determine the quantum fluctuations of 
the signal field in a SPOPO below threshold. Let us 
introduce the following operators: 



m 

S[n,k(t) ^ ^ Lfe m 5i nyn (t) . 



(34) 
(35) 



As Lk ■ Lk> = J2m Lk,mLk' 



3fe,fc' 



one has trivially 



S k (t),Sl(t)\ = 5 kfk >, 

'SinAt),SLk,(t')]^5k,k'5(t-t'), 

and the correlation 

[§*A*)>8LvV)) = 8k,k>S(t-t), (36) 



as well. Hence Sk and Sin,*; are the annihilation opera- 
tors of a combination of signal modes of different frequen- 
cies, which are the eigenmodes of the linearized evolution 
equation (|26[) . The corresponding creation operator ap- 
plied to the vacuum state creates a photon in a single 
supermode which globally describes the frequency comb, 
or train of pulses. Analogously one can define supermode 
output operators 



Sout,k(t) — ^ Lk .mSout.m (t) , 



(37) 



where the output boson operator s ou t, m (t) relates to the 
intracavity and input boson operators through the usual 
input-output relation of high finesse optical cavities, 



^out,m (0 ^in.m 

One can then write: 
d 



(t) + v/2^s m (t) 



(38) 



dt 



-7sSfe + %aA k Sl + \/2<y s S inik . 



(39) 



Let us now define quadrature hermitian operators S k 
by: 



(±) 



S k + — Sk + s| , 



s. 



(-) 



(Sk-sl), 



(40) 
(41) 



and analogously for Si n! fc and S ou t,k, which obey the fol- 
lowing equations: 



dt 



AW 



(42) 



with At given by Eq. (13T))) . These relations enable us 
to determine the intracavity quadrature operators in the 
Fourier domain S^ (ui) 



iwS^H = A^S^M + v^-SgfcH- ( 43 ) 
Finally, the usual input-output relation on the coupling 
mirror (1381). which can be written as 



= -Si n , m (w) + yj2rf s S m (u), 



(44) 



being Sout.m^) the Fourier transform of the output bo- 
son operator s ut,m (t), extends by linearity to any su- 
permode operator as the mirror is assumed to have a 
transmission independent of the mode frequency. One 
then obtains the following expression for the quadrature 
component in Fourier space of any signal supermode, 



S. 



(±) 

out,fc 

4 ±} H 



(±) 



7 s (l±rA fc /A ) -iuj 



k K ' 7 s (-l±rAk/A )+wj' 
One has also, for the operators in Fourier space: 



(45) 
(46) 



(6) 



V 



(o,6) 



2tt 
±,b=± 



S M 6 (ui + w 2 ) (47) 
(48) 



with t/ +:+ ) = tj' ' ' = 1 and rf + ' ) = —rf ,+ ) = i. 
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B. Homodyne detection 

The variances of the quadrature operators can be 
measured using the usual balanced homodyne detection 
scheme: the local oscillator (LO) is in the present case a 
coherent mode-locked multimode field El (t) having the 
same repetition rate as the pump laser: 



E L (t) = E^(t)+E+(t) 
#l (*) = *<! e ™ e ~^'" 

i n 

Et{t)=[EZ{t)]\ 



(49) 
(50) 

(51) 



where J2 m \ e m\ 2 = 1, and £l is the LO field total am- 
plitude factor. The output signal exiting the SPOPO, 
E s ,out (t), is combined with El (t) in a 50%-50% beam 
splitter, the intensity of the two output ports is measured 
using photodiodes of unity quantum efficiency and their 
difference constitutes the homodyne signal. Writing 

-E s ,out (t) = E s out (t) + E+ out (i) , 

^"s.out W = *£out / Sout,m (t) e b,m , 



^s.out (*) - [ E s,out (*)J J 

where £ ut is a proportionality constant. If sufficiently 
fast detectors were used the measurement would give an 
instantaneous signal represented by the operator 

When detectors are not so fast (we are considering inter- 
pulse separations on the order of few ns) they average 
over many pulses along their response time and i must 

be substituted by ?h (t) = ^ It-ra/2 ^'*(*0) wmcn can 
be very well approximated by 



E£ (*) E+ out (t) + E+ (t) E~ out (t) 



*H (*) =J2[ 



"i°out,m 



(52) 



where we considered that Td ^> 2tt/Q and used that 

^out.m 

(i) and §l ut m (t) vary little during the time r d 
[11], what roughly requires that Td <C 7 S _1 . Note that 
this case, namely 2tt/Q -C Td -C 7s" 1 , is sensible as 
7 S 7 1 = 2T.T 1 (2tt/H), where T s -C 1 is the transmission 
factor of the single cavity mirror at which signal losses 
are assumed to be concentrated. Then operator in (jB"2"|) 
represents the outcome of a balanced homodyne detec- 
tion that uses as a local oscillator a modelocked laser with 
the same repetition rate as the SPOPO and with spectral 
components given by e m . The variance of «h measures 
then the fluctuations of the projection of the output field 
on the local oscillator. 



1. Perfect mode matching case 

When the coefficients e m of the LO field spectral de- 
composition are equal, apart from a global phase <f>L, 



to the coefficients Lk t m of the fc-th supermode, e m = 
e*^ L £fe m , one measures, according to (|52"]l . a photocur- 
rent difference proportional to 

Jh (*)=e^sL > fc(*)+e- < * L Sout,*(t) 

= O) cos^l + Sin ^- (53) 

The two following variances, depending on the local 
oscillator phase value 0l, are measured (see next subsec- 
tion for the demonstration), 



,-<->, , {-) { \ (.-)/ \ 7s (l-rA fc /A ) 2 +cj 2 
^ M = ^ H «i (-w) = — - 2 , 

7^(l + rA fe /A ) 

(54) 



,(+) r,,^,(+) 



(55) 



where vj^ are given in (|46[) . Equations (|54I55[) show 
that the device produces, as expected, a minimum un- 
certainty state and that quantum noise reduction below 
the standard quantum limit (equal here to 1) is achieved 
for any supermode characterized by a non-zero A^ value. 
Clearly which quadrature is squeezed depends on the sign 
of Afc/Ao, so that when positive, it is Sf. ' the squeezed 
quadrature (phase-quadrature squeezing) and vice-versa 
(amplitude-quadrature squeezing) . The smallest fluctua- 
tions are obtained close to threshold (r = 1) and at zero 
Fourier frequency (oj = 0): 



(Vfc) • = 



A - |A fc 
A + |A fe 



(56) 



In particular, if one uses as the local oscillator a copy 
of the critical mode k = (identical to the one oscil- 
lating just above the threshold r = 1) one then gets 
perfect squeezing just below threshold and at zero noise 
frequency, just like in the c.w. single mode case. But 
modes of k ^ may be also significantly squeezed, pro- 
vided that |Afc/Ao| is not much different from 1. In the 
next Section we analyze the behavior of the squeezing 
levels just described. 

Our multi-mode approach of the problem has therefore 
allowed us to extract from all the possible linear combi- 
nations of signal modes the ones in which the quantum 
properties are concentrated. 



2. General case 

As always in quantum optics, the measurement of a 
high degree of squeezing in SPOPOs requires the use 
of a mode matched LO, namely of spectral components 
e m = e^ L Lfe im . It is not always an easy task and was 
recognized in [l[ as the main experimental limitation in 
pulsed squeezing. With the present ultrashort pulses, 
one can use pulse shaping techniques with the help of 
dispersive elements and programmable phase modulators 
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39, 40, 41]. In view of future experiments, it is important 
to determine the noise levels measured using a LO of ar- 
bitrary shape, in order to know the accuracy with which 
the perfectly modematched LO must be approached us- 
ing pulse shaping techniques. We derive in this section 
the noise spectrum for a LO of arbitrary shape, that we 
will use in section VI. 

Let us define the projections of the LO frequency comb 
onto the supermodes as 



m = 



m 

This expression can be inverted to yield 



(57) 
(58) 

where the well known result Lk im Lk, n = &m,n involv- 
ing the elements of a basis has been used. Substitution 
of Eq. jSHD into Eq. O) yields 



z H (t) = ]T [ Re + Im ( d *0 sLUt)} , (59) 

fc 

where the quadrature operators, Eq. I|40p. have been 
used. 

The noise variance spectrum associated to zh (t), 
V (w), can be computed as 

/+oo 
dr(J H Wfe(< + r))e- !wr . (60) 
-OO 

Using Eqs. (gSJ), flU, |g7J) and (J52J, it is finally equal 



to: 



+i Redfe Imdi 



,(-) 



m4 +) n]} 



(61) 



Equation ([61]) gives the general expression of the squeez- 
ing spectrum corresponding to a generic LO defined by 
its supermodal amplitudes dk given by Eq. (|57p. When 
the LO is proportional to the supermode labeled by k, 
say e m = e 4 ^ L ifc,m, and </>l = 0,7r/2, the two special 
quadratures (|4"tJ)) arc selected and the results (|54"|) and 
(l55l) are recovered. 



V. DIAGONALIZATION OF THE MATRIX C: 
ANALYTICAL AND NUMERICAL RESULTS 

We have seen that all the properties of the SPOPO are 
directly related to the series of eigenvalues A*, depending 
both on the phase matching properties of the crystal and 




FIG. 1: (Color online) Sketch of the phase-matching matrix 
,fm,q in the space of the integer numbers {m, q} corresponding 
to frequencies {uj m ,Lu q }. In white are represented the regions 
of optimal phase-matching while in blue the unmatched re- 
gions. 



pump spectral characteristics. We will consider now in 
more detail the characteristics of these eigenvalues. 

Fig. Q] shows schematically the typical appearance of 
the phase-mismatch matrix f m _ q — sl ^ m - q (Eq. ([22])). 
for a typical configuration (see Section VI for details and 
real examples). It has the shape of a hyperbola, whose 
branches, of width Ni, display a minimum distance be- 
tween them called d. Another relevant quantity is the 
"width" N2 marked in the figure. These quantities will 
be useful for determining the Gaussian limit, that we will 
consider in the next section. 

The matrix C is the product of f m ,q with the pump 
spectrum a m+q . As the last quantity is constant for 
to + q = constant, the pump selects a portion of ma- 
trix f m .q, roughly given by the intersection of f m<q with 
a straight band oriented along the direction m + q = 
2TO max , where TO max corresponds to the maximum of the 
pump spectrum. 



A. Analytical approach 

When the pump spectrum is not very broad, the re- 
sulting nonzero matrix elements of C are confined within 
an "ellipse" whose principal axes are oriented along the 
directions to + q = and to — q = 0. In this case one can 
forget the secondary maxima of the sine function and use 
a Gaussian approximation for the phase-matching ma- 
trix: 



lit 



(62) 



The form for f m q follows from the approximations ^is^ 
e~~ and mn £ 2 ~ e 



e ^2 [27J, where the parameters 771 
and r]2 can be opportunely chosen so that the results 
from the diagonalization of the coupling matrix obtained 
from (|62|) match optimally to the results of the numerical 
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diagonalization of £ m ,q- By choosing them as r\\ = 5 and 
r/ 2 = 12, the following expressions for the widths Ni and 
N2 can be obtained: 



No = 



2\/3 



(63) 
(64) 



where the (5 coefficients are the ones introduced in (|23l) . 

iU, (EHD. 

Let us assume in addition that the pump has a Gaus- 
sian spectrum centered at 2lu : 



where \a m \ 2 — 1 and iV p 
the number of pump modes, 
are then equal to: 



2 e 2 V«pJ (65) 

= (fiT p ) 1 is a measure of 
The matrix elements C m ,q 



(note that we changed the notation from indices to argu- 
ments) . 

The eigenvalues are given by 



A fe = A p k , 



where 



A = ir 1/4 ^2N~ p 



■1 + 2, 



and 



n = 



|^p 



10 



T2 



4V3 



(70) 

(71) 
(72) 



(73) 



N 2 



(66) 

The Gaussian approximation (|66p is correct as far as 
the interplay between the pump spectrum and the phase- 
matching is opportune. More quantitatively speaking, we 
have to demand that the two branches of the hyperbola 
in Fig. [T] are sufficiently separated each other, which 
corresponds to require that d ^> Ni (we will consider d > 
lOiVi), and the pump width is sufficiently smaller than 
the width of phase-matching function along the direction 
m — q = 0, which corresponds to the condition N p < Ni. 
These conditions lead to the following bounds: 



PI 



I /?2p — \@2s 



> 20, 2 |/3 X 



N p < 1. 



In terms of the crystal parameters and of the pump pulse 
duration these validity limits can be cast as 



for. 

Z > 20 1 p 



1 h"\ 

2 K s I 



(k' p fcj 
> \k'-k'\l. 



(67) 
(68) 



For a BIBO crystal under typical conditions the above 
inequalities read I > 0.2 /im (hence it is not a serious 
condition) and r p > 400 fs x I /mm. Hence, for a crystal 
length Z = 1 mm, t p should be larger than 400 fs in order 
that the Gaussian approximation is valid, while for I 
0. 1 mm the condition is met just for r p > 40 fs. We note 
that condition (|67[) means that the pump duration should 
be longer than the temporal walk-off between signal and 
pump modes along their propagation inside the crystal. 

The eigenvalues of such a matrix turn out to have a 
simple analytical expression at the continuous limit, i.e. 
when one can replace in (|29[) the sum by an integral, so 
that: 



A k L k (m) 



dqC (m, q) L k (q) 



(69) 



Ao and p are given in the limit ti <SC t\, which holds 
unless the crystal length I < 0.1//m; hence p is very close 
to —1. Equation (|70|) corresponds then to an alternating 
geometric progression of ratio p, whose first element Ao 
is positive. 

The eigenvectors are similar to the well-known 
Hermite-Gauss TEM pg transverse modes. They are given 
by 



2\nJ jj. [ 



(74) 



where H k is the Hermite polynomial of order k, and A^ s 
is the number of signal modes. Hermite-Gauss functions 
being simply proportional to their Fourier transforms, 
their temporal shape is exactly the same as their spectral 
shape. The pulse duration of the zeroth mode, r s , is given 
by 



= 2T 2 J T f+T, 



P ' 



(75) 



Under typical conditions [42| the times t\ and T2 are on 
the order of t\ ~ 100 fs and t 2 ~ 5 fs for a crystal length 
I = 1 mm, and in general t 2 *C t\ whenever I > 0.1 pxa.. 
Note that the condition (|67p implies that r p 3> rf, so 
that Al ~ 2y/^N p . 

We are then led to the important conclusion that, ac- 
cording to Eq. (|3"Tj) . the SPOPO threshold is roughly 
equal to the cw single mode threshold Pq divided by the 
number of pump modes, and can be therefore very low. 
For example, if N p = 2 x 10 4 (corresponding to r p = 100 
fs and a cavity length L = 2 m) and considering the 
case already discussed (a 100 /im-thick BIBO based lin- 
ear SPOPO pumped at 0.4 pm), we expect, for T s = 0.01, 

a pump irradiance at threshold of 0.1kWcm~ 2 , 

and an average pump power of 16 mW for a typical pump 
beam radius of 70 pm. 
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B. Numerical approach 

In the general case one must diagonalize numerically 
the 10 5 x 10 5 matrix C. The situation can be dramatically 
simplified from the computational viewpoint by noting 
that a scale transformation affecting the SPOPO param- 
eters allows diagonalization of a much smaller matrix. 

Let us now consider a set of parameters defined by 

(3[ = k/3i, (3' 2p = k 2 (3 2p , /4 = k 2 (3 2s , (76) 

N p = n^Np (77) 

with k a large and positive real number. Let us call 
£' (m, q) the value of the matrix element with these new 
parameters. The form of the matrix coefficients C (m, q) 
when the phase mismatch coefficient <j) m ,q has been re- 
placed by its approximate value (|22[) implies that : 

C! (m, q) = \[kC (nm, nq) . (78) 

Let us set the eigenvalue problem for C': 

A' k L' k (m) = J dqC (m, q) L' k (q) , (79) 

where we added a prime to denote the new eigen- 
elements. Using (1781 and performing the change of vari- 
ables x = nm, y = nq, one finds that: 

A fc = V^K, (80) 
L k , m = L' k (m//s) . (81) 

These two relations are very useful as they allow to com- 
pute numerically eigenvalues and eigenvectors of C in 
terms of the corresponding ones of much smaller matrix 
C because, according to (|78[). the support of C is much 
reduced as compared with that of C. In any case the 
value for k must be chosen adequately in the sense that 
the diagonalization of the toy problem can be cast in the 
integral form (|T9"|) so as to keep £ a smooth function of 
(m,q). 

As a by-product of the demonstration an interesting 
prediction on the influence of the cavity length can be 
drawn: consider that, given a SPOPO, we modify its 
length according to L' = kT x L. This modifies the free 
spectral range as ft' — k£1 and the new SPOPO param- 
eters relate to the old ones as in ([76]) and (f77[) . Hence 
(1791 leads to 

A' fe = v/ZyZAfc. (82) 

This is the case in particular for Ao and the new pump 
threshold becomes P t ' hr = (L/L')P t h r - Hence increas- 
ing the cavity length (and correspondingly decreasing the 
repetition rate) decreases the threshold accordingly. 

VI. APPLICATION OF RESULTS IN 
REALISTIC CASES 

In this Section we discuss the threshold and squeezing 
properties of experimentally realizable SPOPOs. This 



study requires the numerical diagonalization of the ma- 
trix C. We have explored many different configurations 
involving different pump pulse durations r p , different cav- 
ity lengths L, different crystal thicknesses I and even dif- 
ferent phase-matching conditions (critical and noncriti- 
cal) that give rise to different dispersion properties. We 
have considered both BIBO and KNb03 crystals and 
have obtained similar results in the sense that the an- 
alytical approach given above describes very well what is 
numerically found in the region (|67p . no matter the par- 
ticular values of the parameters. When that condition 
gets violated, deviations from the analytical result are 
obviously found but they affect mostly the behavior of 
the eigenvectors, not so much the one of the eigenvalues. 
As the analytical limit is the best also from an experimen- 
tal viewpoint (there the eigenvectors are Hermite-Gauss 
modes, which can be reasonably easily produced exper- 
imentally) we consider here one case that clearly fulfills 
condition (|67[) with parameters compatible with the ex- 
perimental setup that is currently under preparation. For 
the sake of completeness we also consider another one 
that "slightly" violates condition (f6"7[) . We wish to re- 
mark that these cases are representative of what we have 
found in an exhaustive study. Finally, when condition 
(f6"7j) is more severely violated large deviations from the 
Hermite-Gauss case are observed that give rise in fact to 
new phenomena that deserve a study on their own and 
are not treated here. 

The cases we discuss here correspond to collinear, de- 
generate type I critical phase-matching at 0.4/zm pump- 
ing of a BIBO crystal [33j |. obtained when the pump po- 
larization is ordinary (parallel to the direction Ox) and 
that of the signal is extraordinary (o — > e + e). Us- 
ing Sellmeier's coefficients for BIBO we obtain that such 
phase-matching occurs at an angle 8 — 151° between 
the direction Oy and the direction of propagation of the 
pump (and the signal) beam, in agreement with |43| . For 
this configuration we obtain the values for the dispersion 
parameters reported in Table [Til Also given in that ta- 
ble are the values of the free spectral range f2 and pump 
pulse duration r p that will be used along this section. We 
shall assume a pump with Gaussian spectrum and cen- 
tered at the phase-matched frequency 2ujq. Finally we 
consider three different values of the crystal length: 

case A: I — 0.1 mm, 
case B: / = 0.5 mm, 
case C: I = 5 mm, 

With all these values one can compute the phase mis- 
match angle <p in the second order dispersion approxima- 
tion, see Eq. (|2"2"]) , and finally the matrix C. 

Case A verifies well the condition l|67p . which reads 
here r p > 40 fs. On the contrary cases B and C do not 
verify it, which now reads r p > 200 fs and t p > 2000 fs, 
respectively. 



10 



(a)BIBO dispersion parameters 



(b)Pump 
parameters 





k' (s m- 1 ) 


k" (s 2 m- 1 ) 


fi (MHz) 


T P (fs) 


pump 


6.6537 x 10~ y 


4.7248 x KT 2b 


2tt x 75 


100 


signal 


6.2664 x 10" a 


1.6420 x 10" 2b 



TABLE II: Dispersion parameters for BIBO crystal and 
pumping field. 



A. Case A 

In the next figure we show the matrix C corresponding 
to this case. Fig. [2] corresponds to the frequency (inte- 
ger indexes) representation, which is the actual matrix 
to be diagonalized. Phase matching occurs in the lighter 
regions. Darker regions are highly phase-mismatched. 
Results of the numerical diagonalization, obtained as ex- 
plained in Section l"VB I by using a scale factor k — 1000, 
are shown in Figs. ([3]) and (j4|). 

In the singly resonant threshold of Wthr — 

29 mW is readily obtained from Eq. (|3l"|) for the cor- 
responding eigenvalue Ao ~ 270, by considering a trans- 
mission factor of T s — 0.01 and a beam waist of 70 /im. 
This results are in perfect agreement with the analytical 




FIG. 2: (Color online) Case A: C matrix. 



< 




FIG. 3: (Color online) Case A. Spectrum of eigenvalues. 
Comparison between the numerical and analytical solutions. 
Among about 10 5 supermodes only a relatively small part 
(~ 40) is dynamically significative. 




0.05p 
0.025 
- 



m (x10 b ) 

(b) 



0.05 
0.025 



-10 -5 5 10 

m (x10 5 ) 




(c) 



FIG. 4: (Color online) Case A. Comparison between the nu- 
merical (solid blue line) and the analytical (dotted red line) 
solutions of the eigenvectors associated to the four highest 
|A*|. 



value predictable by the expression Eq. (|7lj) . For a dou- 
bly resonant cavity this result has to be multiplied by 
means the correction factor that accounts for the geom- 
etry and is reported in Table |TJ 

From Eqs. and we can calculate also the 

noise reduction corresponding to the first four eigen- 
vectors shown in Fig. 4 for a zero noise frequency. 
Evidently we are assuming to be able to master the 
spectral shape of the local oscillator in order to ex- 
actly match it to the supermode whose noise variance 
spectrum is to be measured. Nevertheless, the opti- 
mization of mode-matching between the local oscilla- 
tor and a specific supermode can result a difficult task 
even when pulse shaping techniques are used. Let's con- 
sider, then, the case where the best we can do is to deal 
with a local oscillator shaped as a Gauss-Hermite poly- 
nomial e Km = 7 r- 1 / 4 AT i ; 1/2 e-5( m /^L) 2 e ^H k (m/N L ), 
where JVl is the number of longitudinal modes of the lo- 
cal oscillator comb. In such situation the variances have 
to be evaluated using the general expression Eq. (|6l"j) . 
where the noise variance spectrum is given by the sum of 
all the supermodes noise variance spectra weighted by the 
mode matching parameters d^, which describe how well 
each supermode projects on the local oscillator field. In 
Table lTlTl we compare the degree of squeezing measured in 
the situation of perfect mode-matching and the situation 
where the best local oscillator is Gauss-Hermite function 
of adjustable spectral width, for a pumping power 20% 
below threshold (i.e. r = 0.9). In the latter case, the 
minimum noise is obtained around — 2.2 x 10 5 . Such 
comparison evidences the fact that the differences be- 
tween the two situations are small. Hence, the exact 
knowledge of the supermodes shape is not necessary and 
a good degree of mode matching can be obtained simply 
controlling the spectral width of a Gauss-Hermite local 
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V k (dB) 


k = 


k = 1 


k = 2 


k = 3 


perfect 


-25.6 


-24.9 


-24.1 


-23.3 


G-H 


-25.5 


-24.8 


-23.8 


-22.6 



TABLE III: Case A. Comparison between the noise variances 
evaluated, at u — and r — 0.9, in the case of perfect mode 
matching of the LO with the supermodes corresponding to 
k = 0, 1, 2, 3 and the case where the LO is the Gauss-Hermite 
(G-H) function described by the spectral amplitudes ek, m and 
A L = 2.2 x 10 s . 

oscillator. Evidently this circumstance is verified as far 
as the condition for the Gaussian approximation of the 
coupling matrix is respected. However, the number of su- 
permodes that present marked quantum characteristics 
results to be greater than four. By considering, in a qual- 
itative way, that — 5dB is a still significative degree of 
squeezing, we found that all the supermodes correspond- 
ing to the first 45 higher values of |Afe| have variances 
smaller than the considered bound. This is an impor- 
tant result since it proves that SPOPOs are multi-mode 
sources of non-classical light. 

B. Case B 

In the next figure we show the matrix C corresponding 
to this case. Fig. [5] corresponds to the frequency (integer 
indexes) representation, which is the actual matrix to be 
diagonalized. Phase matching occurs in the lighter re- 
gions, while darker regions are highly phase-mismatched. 

m (x10 5 ) 




-250 ' 1 1 1 1 1 

25 50 75 100 125 150 

k 

FIG. 6: (Color online) Case B. Spectrum of eigenvalues. 
Comparison between the numerical and analytical solutions. 
Among about 10 s supermodes only a relatively small part 
(~ 23) is dynamically significative. 




m(x10 5 ) m(x10 5 ) 



(c) (d) 




FIG. 7: (Color online) Case B. Comparison between the nu- 
merical (solid blue line) and the analytical (dotted red line) 
solutions of the eigenvectors associated to the four highest 
|A*|. 



FIG. 5: (Color online) Case B: £ matrix. 

Results of the numerical diagonalization, obtained as 
explained in Section lVBl bv using a scale factor k = 1000, 
are shown in Figs. [6] and [7] 

In the singly resonant threshold of Wthr — 

1.5 mW is readily obtained from Eq. (f3l"j) for the cor- 
responding eigenvalue Ao ~ 235, by considering a trans- 
mission of T s = 0.01 and a beam waist of 70 fim. The 
threshold obtained from the analytical solution is about 

1.6 mW, which is not too much different from the exact 
solution. Hence, even if the experimental situation con- 
sidered here does not strictly verify the conditions for 
Gaussian approximation, we find still a good agreement 



between the numerical and the analytical predictions. A 
qualitative statement about the good agreement between 
the numerical and analytical solutions can be grounded 
also from the comparison between the eigenvalues shown 
in Fig. El 

Assuming perfect mode matching, from Eqs. (|54p and 
(153)1 we can calculate also the noise reduction correspond- 
ing to the first four eigenvectors shown in Fig. [7] at 
the carrying frequency and 20% below threshold (i.e. 
r = 0.9) and compare it to the general case where the lo- 
cal oscillator is described by the Gauss-Hermite spectral 
amplitudes ek, m - In this case, the detection is optimized 
for a spectral width of about Ni, ~ 1.2 x 10 5 . The results 
are reported in Table HVl In this case, the amount of 
squeezing detected using a Gauss-Hermite local oscilla- 
tor, despite of the optimization of its spectral width, is 
not as much as the perfect case, even if still significative. 
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V k (dB) 


k = 


k = 1 


k = 2 


= 3 


perfect 


-25.6 


-24.1 


-22.6 


-21.0 


G-H 


-23.9 


-20.8 


-17.9 


-15.0 



TABLE IV: Case B. Comparison between the noise variances 
evaluated, at u — and r — 0.9, in the case of perfect mode 
matching of the LO with the supermodes corresponding to 
k = 0, 1, 2, 3 and the case where the LO is the Gauss-Hermite 
(G-H) function described by the spectral amplitudes ek, m and 
N L = 1.2 x 10 s . 

This result is a consequence of the fact that the situa- 
tion now considered is slightly violating the bounds for 
the Gaussian approximation and, hence, the supermodes 
have spectral amplitudes that are no more characterized 
by Gauss-Hermite functions. Nevertheless, even if Gaus- 
sian approximation is not perfect, its prediction capabil- 
ity is still relevant. 

Since we are interested to SOPOs as multi-mode 
sources for non-classical light, let's consider the same 
qualitative argument we considered in the previous sec- 
tion. In this case about 23 supermodes present a degree 
of squeezing greater -5 dB. Despite the fact that, with 
respect to the O.lmm-thick crystal, the number of super- 
modes that characterize the SPOPOs output is smaller, 
it can still be considered highly multi-mode. 

C. Case C 

The last case we consider corresponds to a configu- 
ration that is strongly non-Gaussian. This can be di- 
rectly observed by comparing the matrix C obtained in 
this case and reported in Fig. [8] with the matrices for 
cases A and B. In fact, since the pump pulse duration 



m (xio 5 ) 




FIG. 8: (Color online) Case C: £ matrix. The inset is the 
magnification of the matrix around the phase-matched fre- 
quency corresponding to m = 0, q — 0. 

t p is much smaller than the temporal walk-off \k' — k' s \l, 
in the space {m, q}, the pump bandwidth N p is larger 
than the phase-matching bandwidth Ni (see Fig. [!}. 
Consequently, the pump selects not only the principal 
peak of the "sine" function corresponding to the phase- 
matching matrix f m ,q, but also several secondary max- 
ima. This has an important consequence for what con- 




50 100 150 200 

k 



FIG. 9: (Color online) Case B. Spectrum of eigenvalues. 
Comparison between the numerical and analytical solutions. 
Among about 10 s supermodes only a relatively small part 
(~ 125) is dynamically significative. 




(c) (d) 

FIG. 10: (Color online) Case C. Comparison between the 
numerical (solid blue line) and the analytical (dotted red line) 
solutions of the eigenvectors associated to the four highest 
|A*|. 



cerns the eigenvectors and eigenvalues of C, obtained as 
explained in Section IVBl by using a scale factor k = 200, 
that we report in Figs. l9land[l"0l 

In Fig. O the spectrum of the eigenvalues obtained 
from the analytical solution (|T0[) (red circles) shows a 
great discrepancy with the eigenvalues obtained from nu- 
merical diagonalization of C, as expected. In particular, 
the part of spectrum, corresponding about to the first 
50 eigenvalues, flatten around the critical value Ao ~ 36, 
while the Gaussian approximation predicts always a geo- 
metric progression-like behavior with a critical eigenvalue 
A Gau SS _ 44 In the 

same experimental configuration 
as previous cases (singly resonant cavity, transmission at 
signal frequencies of T s = 0.01 and beam waist of 70 (Jxa), 
a threshold of Wthr — 0.67 mW can be obtained from Eq. 

(ED- 
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V k (dB) 


k = 


k = 1 


k = 2 


k = 3 


perfect 


-25.58 


-25.58 


-25.57 


-25.57 


G-H 


-25.27 


-24.28 


-22.51 


-20.22 



TABLE V: Case C. Comparison between the noise variances 
evaluated, at u — and r — 0.9, in the case of perfect mode 
matching of the LO with the supermodes corresponding to 
k = 0, 1, 2, 3 and the case where the LO is the Gauss-Hermite 
(G-H) function described by the spectral amplitudes ek, m and 
At = 0.08 x 10 5 . 



The result of a threshold higher than the one expected 
in the Gaussian approximation has a physical explana- 
tion. Since, for the time durations involved, the process 
of parametric down conversion can be considered almost 
instantaneous, for a mode-locked pumping field the peak 
power necessary to reach the oscillation threshold is the 
result of the coherent contribution of all its modes. This 
circumstance is formally expressed by the fact that the 
analytical expression for Ao (see Eq. (17(H) ) depends on 
the number of modes in the pump pulse N p . But, beyond 
the Gaussian limit, the fact that N p ^> Ni implies that 
not all the N p pump modes are equally phase-matched 
and, then, not all can optimally transfer energy towards 
the signal modes. The same phenomenon can be under- 
stood even in the temporal domain. In fact the quantity 
|fcp — k' s \l corresponds to the temporal walk-off accumu- 
lated by the pump and signal pulses through a passage in 
the nonlinear crystal. When, in a non-Gaussian regime, 
the condition ([6"Tj) is violated, the walk-off between pump 
and signal is bigger than the pump width r p and the two 
field cannot optimally exchange energy all along the crys- 
tal length thus increasing the instantaneous peak power 
necessary to reach the oscillation. 

On the other hand, in Fig. [TO] the eigenvectors re- 
trieved in the Gaussian approximation (|T4[) do no more 
fit the numerical solutions, as expected. In fact, even 
if they still preserve a shape similar to Gauss-Hermite 
functions, they result to be shorter in the domain of fre- 
quencies and are affected by a small modulation of the 
spectral amplitude. 

As the previous cases, we consider the quantum prop- 
erties of the supermodes and compare the situation of 
perfect mode matching of the LO with that of a Gauss- 
Hermite LO. In this latter case the detection is optimized 
for a spectral width of about N L ^ 0.08 x 10 5 . The re- 
sults are reported in Table [V] for the eigenvectors corre- 
sponding to the first four biggest |Afc| . The fact that 
they are close to degeneracy (see Fig. \§§ is reflected 
in an almost equal reduction of noise variances below 
the standard quantum limit. Despite the differences re- 
ported between the variances evaluated both by means 
of a perfectly mode matched and a Gauss-Hermite LO, 
a still significative degree of squeezing can be detected 
in realistic situations thus suggesting that the shaping of 
the LO is not a critical issue for the experimental con- 
figuration considered in this section. Actually, there is 
a larger set of supermodes the variances of which are all 



close to the value of the critical one. In particular, there 
are about 30 supermodes that have variances comprised 
in 1 dB, between —24.6 dB and —25.6 dB. Furthermore, 
by considering the number of supermodes that present a 
noise reduction bigger than —5 dB, one discovers that, 
this time, their number amounts to about 125. 



D. Discussion on the influence of the crystal length 

We have seen that even if the cases A and B do not 
verify at the same time the condition (f6"T|) . the analytical 
solution obtained in the Gaussian approximation works 
quite well in both cases. Nevertheless, in spite of the fact 
that this condition gives an approximately good idea of 
the reliability of Gaussian approximation, it is interesting 
to study the passage from a perfectly Gaussian case to 
a non-Gaussian one trough a "gray" region where the 
differences between the two cases are not big. 

Let us consider Eq. (|TD|) in the limit of very large 
/. In such case, since from Eq. (|73[) t\ ^ t p , then Ao 
asymptotically converges to: 

Ao^y^_W (83) 

This expression indicates that the product Aq x I is con- 
stant for values of / compatible with a non-Gaussian 
regime. In Figure 111! then, we report the values of 
this product as a function of the crystal thickness. For 
I < 1mm the analytical solution, as expected, is in 
good agreement with the numerical one, while for greater 
thicknesses the discrepancy is significative. Also this re- 
sult is expected, since the analytical solution for the crit- 
ical eigenvalue has not validity when the condition ([6"T[) 
is violated. On the other side, the fact that also the ana- 
lytical solution reaches asymptotically, for increasing I, a 
plateau suggests a 1 /Z-like behavior of A . The existence 
of such a plateau can be explained from the point of view 
of the evolution of the pump and signal pulses in the time 
domain. As discussed in the previous section, when the 
condition (|67jl is violated, the walk-off between pump and 
signal is bigger than the pump temporal width. As a con- 
sequence, the exchange of energy between the two fields 
is disadvantaged till a point where the threshold can- 
not change anymore even increasing the crystal length. 
Since, from Eqs. (gHJ) and ([31]), P thl cx (A x Z)~ 2 , for 
large values of I, then, also the product Ao x I reach a 
constant value, thus explaining the plateau in Fig. 1111 

In the same way, one can explain the discrepancy ob- 
served between the two plateaux in Fig. [11] As already 
discussed in case C, in non-Gaussian configurations the 
pump bandwidth is larger than the phase-matching one, 
thus not all the pump modes are phase matched and not 
all of them contribute to the final value of the threshold. 
Therefore the number of pump modes actually involved 
is smaller than the nominal value N p that should, then, 
be corrected. 
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FIG. 11: (Color online) Numerical and analytical curves for 
the product Aq x ( versus the crystal length. 



From a quantum point of view, we have detailed in 
the previous sections the noise properties of the super- 
modes connected to the first four highest Aj, (see Tables 
Mil IIVI and [Vj and we calculated the number of super- 
modes showing a squeezing better than —5 dB for getting 
a qualitative indication about the "multimodicity" of the 
system prepared in a specific experimental configuration. 
These results can be appreciated in Fig. [T3] where the 
noise variances of the supermodes that satisfy this cri- 
terium have been traced for the three cases previously 
discussed. The curve in the middle corresponds to case 
A (I — 0.1mm), a Gaussian configuration. As the thick- 
ness of the crystal is increased to I = 0.5 mm (case B) 
the condition (|67[) is violated but the Gaussian approxi- 
mation is still good. This means that, even if the value 
of Ao is decreasing (because we are reducing the number 
of pump modes that are phase-matched) the spectrum is 
still a geometric progression but with a smaller, in abso- 
lute value, common ratio (see p in Eq. (|70p ) thus causing 
an overall decrease of the spectrum with respect to the 
case A. The consequence is a reduction of the number 
of supermodes with a squeezing greater than —5 dB as 
it results from the upper curve in Fig. [T^] (red circles). 
Finally, when the crystal length is further increased to 
I = 5mm (case C), we pass to a completely non-Gaussian 
configuration where the decrease of the critical eigenvalue 
Ao (from ~ 271 for I — 0.1 mm to ~ 35 for I = 5 mm) 
causes a significative deformation of the spectrum and 
a non null set of eigenvalues flatten around Ao. In this 
case, since the degree of squeezing for each supermode de- 
pends on the ratio |A^/Ao| (see Eq. (|5rj|) ). the amount of 
squeezing is globally increased as the lower curve (green 
diamonds) in Fig. H"2l shows. 

These results not only confirm that a SPOPO is a 
highly multi-mode device but also show another impor- 
tant quality: the malleability for controlling its "mul- 
timodicity". We have seen, in fact, that one can just 
increase the thickness of the nonlinear crystal in order to 
improve the number of supermodes that play an impor- 
tant role from a quantum point of view. 



FIG. 12: (Color online) Comparison between cases A (I = 

0. 1mm, middle curve), B (/ = 0.5mm, upper curve) and C 
(I — 5mm, lower curve). Noise variancese evaluated, at u) = 
and r — 0.9, in the case of perfect mode matching of the LO 
with the supermodes corresponding to the significative part of 
eigenspectrum < k < 150. Among about 10 supermodes 
only a relatively small part is significative from a quantum 
point of view. In particular, there are 45, 23 and 125 super- 
modes, respectively, that are squeezed more than —5 dB. 

In this paper we have presented a study of the eigen- 
values and eigenvectors of the matrix £ in function of 
the crystal length I. However, since both the pump pulse 
width and the crystal length are present in Eq. (|67p , no- 
tice that effects similar to those discussed in this section 
can be observed by playing with r p and keeping constant 

1. In fact the passage from a Gaussian to a non-Gaussian 
configuration takes place when the pump bandwidth N p 
becomes smaller than the phase-matching one N\ and, 
clearly, this can be obtained keeping fixed the latter and 
increasing r p (see Fig. [TJ. Eventually one could even 
fix I and r p and play with the group velocity mismatch 
by choosing different types of nonlinearities, but, experi- 
mentally, this could result in a stiffer malleability of the 
device. 



VII. CONCLUSION 

In this paper, we have shown that both the dynami- 
cal and quantum noise properties of the SPOPO depend 
on the spectrum of eigenvalues A& of the linear coupling 
matrix C. We have studied in detail this spectrum in 
various experimentally feasible configurations, using ei- 
ther an analytical approach in some simple limit cases, 
or a numerical approach in the general case. It turns out 
that among the roughly 100,000 eigenvalues, 99,900 or 
so are zero, and about a hundred are significantly dif- 
ferent from zero and contribute to the quantum dynam- 
ics of the system. SPOPOs are therefore devices which 
produce simultaneously many highly squeezed vacuum 
modes. This property can be used to improve the perfor- 
mances of metrological methods using frequency combs, 
for example to perform ultra-accurate time transfer be- 
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tween remote clocks beyond the shot noise limit [46j. In 
addition, it is well known that if one mixes by one way 
or another different squeezed modes, one gets strongly 
entangled states (47| . This is also the case here: we will 
show in a forthcoming publication that SPOPOs are in- 
deed likely to generate various pairs of strongly entangled 
supermodes, as well as multipartite entangled states. 

More generally, this paper is an example of the fact 
that, by using appropriately chosen pump spectra and 
phase matching curves, one can reach various eigenvalue 
spectra, and therefore tailor at will the quantum proper- 
ties of the light generated by the optical system, which 
can be useful for example to generate interesting states 
for multidimensional quantum information processing. 



VIII. APPENDIX : QUANTUM MODEL FOR A 
SINGLY RESONANT SPOPO 



We detail in this appendix the derivation of the Heisen- 
berg equations in the case of a singly resonant degenerate 
type I SPOPO, i.e. when the signal field resonates inside 
the cavity, but not the pump field. We consider here the 
case of the linear Fabry-Perot cavity. The treatment of 
this case is more complex than the doubly resonant case 
usually considered in theoretical approaches as the pump 
field cannot be quantized inside the cavity. 

The signal field inside the nonlinear crystal, which ex- 
tends from z — — 1/2 to z = +1/2, is written as 



E s (z, t) — i ^2 £s, g sin (k s , q z') s q (t) e 



■H.c., (84) 



where z' = z + L/2, and £ s , 9 



Or the 



■il ~ Y e n(u) e , q )A s L ' 

contrary, the pump field is not affected by the cavity and 
hence it is given by a continuum of modes. We shall use 
the common approach of quantizing the pump field in a 
line of length L p with periodic boundary conditions and, 
in the end of the calculations, we will make L p — > oo. We 
thus write 

E p (z,t) = iJ2£(vm)e- iVmt x 

m 

x [pW(f)e^+tI(t)e" M ] + 

+ H.c, (85) 



where the superscripts (±) label the propagation direc- 
tion. £ (V m ) = \J 2e ApLpn J[^j &re sin g le photon field 

amplitudes. In the limit L p — > oo the frequencies v m 
are given by u m — m-jf^, m € N, and the wavenumbers 

K m = n { Vm ' Vm | with n(v rn ) the crystal refractive index. 
Note that we are writing the fields as a superposition 
of plane waves, but the treatment is still approximately 
valid for Gaussian beams as far as the thin crystal is 
placed at the (common) beam waist of pump and sig- 
nal, and the Rayleigh lengths are much longer than the 



crystal length I. In such case Af — irwj with W{ the 
corresponding beam radius. 

The interaction Hamiltonian Hi is calculated as usual 
as in I)?)). Inserting the expressions of the second order 
nonlinear electric polarizations, one obtains the following 
form of the interaction Hamiltonian, in which \ is the 
second order nonlinear susceptibility (whose dispersion 
is neglected): 



Hi(t) 



9 J 



X 



x e i( u) s, m +^ s ,q-Vj)t e ~i(K, m +k Siq )L/2 

+ H.c, 

where we defined the phase-mismatch factor 



pj 

7n,q 



sin [(kj k S:Tn &s,(/) 2 



I 

q) 2 



(86) 



(87) 



In (I86j) we dropped highly phase mismatched terms, as 
usual. 

If we introduce new signal boson operators 
s (t) = s (t) e+ ifc =.™ L /2 



the interaction Hamiltonian becomes as (|86|) but without 
the exponential e - l ( k s,™.+k a , q )L/2 _ j n ^ e f vj owm g we w [\\ 

use the new operators but omit the superscript "new" for 
simplicity. 

From the previous expression of the hamiltonian, one 
can derive following Heisenberg equations governing the 
time evolution of the pump and signal operators: 



dpf ] (t) 
dt 



2h 



x s m (t)s q (t)e- i ^ + ^-^\ (89) 



(ii 



9 3 



x e 



(90) 



The integration of the pump equations yields 

m q 

x / dt's m (t l )s q (t')e^ UJs - m+ ^-"- l '^ t ' , (91) 







where pffj c j — p)^' (0) is the source- free part of the pump 



,(±) 

(the field impinging the nonlinear crystal). 
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Using the usual Wigner-Weisskopf approach, valid be- 
cause the nonlinear interaction is assumed to be instan- 
taneous, and using the approximation 



"2 



(92) 



we obtain a value of ' (t) that one can then insert in 
Eq. (SQ): 



pf (*) = P\Zi (*) + 



s(±) 



m,q 



X S m (f) (t) (5 (cj B ,m + W s ,q - Vj) , (93) 

We can now pass to the continuum limit. For that we 
define continuum pump operators in the following way: 



frequencies u PtT = 2u>q + rf2, rgZ. Thus I\ becomes 

K 1 ) ~ j(w 3im +w a _„)t ' 



m,q 



E 



sin [(fc p> . 



l(Wp,r) 



(^p,r fcs.m ^s,g) 



x / d^ 



(99) 



where fe P)T - = fc (w Pjr ) = fc (2wo + rf2). 
We now define new pump operators 

V^e-^tp^}. (t) = f dve-^pti (v) . 

J r 



(100) 



which can be shown to verify the following property: 



PiJn (*l) Cr 2 (*2) > - <^,r 2 <5 (tl - fc) ■ (101) 



s(±) 



Pfree — \/ 2 7 ^ c P£ree,j J 



p A±) 



(94) 



which verify 



Pfee^i)' Pfee(^) 



a(±) 



Z7TC L„—>oo 



(95) 

Transforming sums into integrals, we obtain the fol- 
lowing equation for the signal modes: 



ds m (t) _ Ai 



dt 



■ix 



£q 
Airhc 



V£ £ sUtM^ 



A 2fr c E E E £s,m£s,q£s,n£s,i 

P n r q 

x s\{t)s n (t)s r {t)I%} q , 



(96) 



where 



m,q 



r(2) 



i/ sin [(fc (z/) - fc Sim - fc s , g ) ■ 
n(i/) (k(v) - k s>m - k s<q ) | 



dz/i5 (w. 



(97) 



n{v) 



*n [(M 



(A; (z^) k s m k St q) 2 



sin - fc s ,» - fc s , r ) 



(/c {i>) k s n k s r ) 2 



x e 



i(ti>8,m+«s,g— f)* 



(98) 



In order to calculate the first integral, we write it as a 
sum over frequency intervals of width £1 and centered at 



Imlq now becomes 



m.q 



^E 



Wp,r Sm [(fcp jr fcs.m ^s,ij) 



< K>,r) (fcp, ; 



h — h 1 i 



X e 



(*)+££(*)!■ (102) 



Retaining only slowly varying terms in the evolution, 
and including the losses of the optical cavity at rate j s , 
one finally gets: 



ds m (t) 
dt 



(*) + Pin rn+q (*) 



(t) + ^J2q s s- m . m (t) + 

q 

9 ^ ' ^ ' fm,qfn,m+q—n X 
n g 

x (t) s n (t) s„ i+g _„ (i) , 
where the coupling constant g is given by 



+ 



(103) 



A\ I / UJQ 



3/2 




and Sin im (t) corresponds to the field at signal frequencies 
entering the cavity through the coupling mirror. When 
that input is coherent or vacuum, the case we consider, 
those "in" operators verify the following correlation 



(sin,™ (t) (*')) = S m,m'S (t ~ t') . 



(105) 



and thus behave as p^' r (t) (see Eq (|10ip ). 

Let us now consider the regime below the oscillation 
threshold : the signal modes are almost not excited and 
the double sum in (|103p can be neglected. Also, the pump 
"in" fields can be approximated by their mean values as 



17 



their fluctuation part gives rise to smaller terms, which 
are neglected for the same reasons as before. Hence, we 
have for a unidirectional pumping : 



< Pin,m (*) > = 







< Phtm (*) >= 



n A p P 
- — / 

2fkoo 



(106) 

P being the average power per unit area of the mod- 
elocked pump laser and X) m l a "i| 2 = 1- Eq. p06[) is 
obtained by demanding that the pump field correspond- 
ing to the set (t)\ equals the external pump field 

given by Eq. (fT]) inside the crystal. The linearized equa- 
tions for the SPOPO below threshold finally become 



ds m CO 
dt 



(t) 



+ 7sO- f-> 



m.q^m-^q^q 



(0. 



(107) 



where 



e c 3 n§T s 



(108) 



32 (xluoY 

In conclusion, we have shown that the linearized equa- 



tions (|107p formally coincide with those of a doubly res- 
onant SPOPO (Eq.(fT5]0. the only difference being the 
exact value of Pa. 
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